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Ori and Thorne have discussed the duration and observability (with LISA) of the transition from 
circular, equatorial inspiral to plunge for stellar-mass objects into supermassive (10 5 — 10 8 Mq) Kerr 
black holes. We extend their computation to eccentric Kerr equatorial orbits. Even with orbital pa¬ 
rameters near-exactly determined, we find that there is no universal length for the transition; rather, 
the length of the transition depends sensitively — essentially randomly — on initial conditions. Still, 
Ori and Thorne’s zero-eccentricity results are essentially an upper bound on the length of eccentric 
transitions involving similar bodies (e.g., a fixed). Hence the implications for observations are no 
better: if the massive body is M = 10 s Mg, the captured body has mass m, and the process occurs 
at distance d from LISA, then S/N < (m/10Mo)(lGpc/d) x 0(1), with the precise constant de¬ 
pending on the black hole spin. For low-mass bodies (m < 7 Mq) for which the event rate is at least 
vaguely understood, we expect little chance (probably [much] less than 10%, depending strongly on 
the astrophysical assumptions) of LISA detecting a transition event with S/N > 5 during its run; 
however, even a small infusion of higher-mass bodies or a slight improvement in LISA’s noise curve 
could potentially produce S/N > 5 transition events during LISA’s lifetime. 

PACS numbers: 04.30.Db, 04.80.Nn, 97.60.Lf 


I. INTRODUCTION 


The gravitational waves emitted during inspiral and 
infall of a body (mass m) into a black hole (mass M) 
should reveal detailed information about the orbital ge¬ 
ometry and the hole’s spacetime geometry, thereby pro¬ 
viding high-precision tests of general relativity [Q . While 
the scattering of waves off background curvature implies 
that waves emitted at any time during the inspiral pro¬ 
vide some small measure of even the smallest scale vari¬ 
ations in the background spacetime geometry, the waves 
emitted as a particle passes through a region provide the 
most sensitive tests of that region: they reveal what path 
the particle has followed, and therefore constrain the 
spacetime to permit such a path. Therefore, to provide 
a sensitive probe of the innermost regions of black-hole 
spacetimes, we want to study orbits that pass as near 
as possible to the hole itself. Unfortunately, this means 
that the signals that are potentially the most informa¬ 
tive about the hole’s innermost structure are typically 
the briefest: they arise from the end of the bound por¬ 
tion of the orbit and from the transition from inspiral to 
plunge. Since the relevant fraction of the orbit persists 
for only a small fraction of the overall detectable inspiral, 
we have significantly less probability to resolve waves dur¬ 
ing this interval than to resolve earlier, longer portions of 
the inspiral. One therefore wants to roughly characterize 
the waves emitted during these intervals (in the case of 
LISA sources, the goal of this paper). If this characteri¬ 
zation suggests that planned observatories such as LIGO 
or LISA could detect them, one should then carry out 
much more detailed studies of these last few orbits and 
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the waves they emit. 

For inspirals appropriate to the LIGO band (~ 10-10 3 
Hz) and which LIGO can plausibly detect (77 = nn/M 
~ 0.01 to 1 ), order-of-magnitude computations (say, by 
post-Newtonian methods) suggest the last few waves are 
detectable ||. But because in this regime simple approx¬ 
imation techniques (such as post-Newtonian j^, Q or 
test-particle approximations) break down, and because 
numerical relativity codes remain incapable of evolv¬ 
ing orbits accurately enough to find the waves, the com¬ 
munity does not yet possess a waveform trusted for any 
purpose beyond detection. 

LISA’s band (~ 10 _ 3 -10 _ 1 Hz) will prove more sensi¬ 
tive to extreme mass-ratio infalls — that is, to stellar- 
mass black holes, white dwarfs, and neutron stars falling 
into supermassive \M = O(10 5 — 10 s ), so 77 ~ 10 -4 to 
10 -8 ] black holes [f|. With such extreme mass ratios, the 
computation of detailed waveforms for purposes beyond 
mere detection should prove much simpler: to under¬ 
stand evolution, we need do nothing more than solve the 
classical radiation-reaction problem, albeit on a curved 
spacetime and with a gravitational, rather than electro¬ 
magnetic, field |t|. While this problem hasn’t been solved 
to the accuracy required to construct long-integration- 
time coherent detection templates, one can employ adia¬ 
batic approximations to address most preliminary inves¬ 
tigations. For example, as Ori and Thorne || have dis¬ 
cussed in the context of circular inspiral, to understand 
the 77 <C 1 transition’s duration — measured in experi¬ 
mentally observable gravitational wave cycles — we do 
not need a precise knowledge of the reaction force. An 
averaged reaction force —- one we can easily deduce from 
the radiation of conserved constants —- suffices for the 
short interval we will coherently employ it. Applying this 
reaction force, we can follow the particle through transi¬ 
tion and thereby predict roughly how long this transition 




2 


will last. 

The goal of this paper is to extend the Ori-Thorne 
analysis to eccentric Kerr orbits, in an effort to estimate 
the prospects of LISA detecting a transition from inspiral 
to plunge. 

This analysis relies on using the radiation of two con¬ 
served constants E, L to compute the effect of radiation 
on the orbit. But for Kerr inclined orbits there is an addi¬ 
tional constant — the Carter constant — whose evolution 
has not yet been related to fluxes at infinity. Since we 
lack the necessary tools, we leave the Kerr inclined case 
to a future paper. 


Outline of this paper and summary of 
conclusions 


In Sec. jn| we will outline the basic physical framework 
behind our approach. In particular, we will introduce 
an explicit procedure to estimate the time duration of a 
transition. This procedure takes as input the net (time- 
averaged) fluxes of energy and angular momentum from 
the particle’s instantaneously geodesic orbits, input one 
obtains from a solution of the Teukolsky equation given a 
geodesic orbit as source. This procedure also takes as in¬ 
put some observationally-defined interpretation of what 
“the transition region” is. As the latter is ambiguous, 
and depends on exactly what sorts of templates one uses 
to find it, the exact length of the transition will depend 
on the convention one uses. 

Ideally, one should define some unambiguous set of 
templates and match those against the simulated emit¬ 
ted waves to both define the transition duration and de¬ 
duce the resulting signal-to-noise ratio for a given source. 
But for brevity and simplicity, as discussed in Sec. Ill, 
we will use a much cruder scheme — based on a purely 
sinusoidal, quadrupolar model for the waves — to char¬ 
acterize the expected LISA signal-to-noise ratio from a 
specific transition crossing. Given S/N for an event and 
loosely-understood rates for transition events, we then 
develop, in Sec. Ill, a scheme for estimating the proba¬ 
bility that LISA will see an event with S/N greater than 
some detection threshold. 

With this complete scheme for estimating the signal- 
to-noise of a characteristic source and determining the 
probability that LISA, in its currently-planned configu¬ 
ration, will see something, in Sec. IV and Sec. |v| we will 
apply it to inspirals into Sc hwar zchilcl and Kerr holes, 
respectively. We find in Sec. VD that Ori and Thorne’s 
zero-eccentricity results are essentially an upper bound 
on the length of eccentric transitions involvi ng s imilar 
bodies (e.g., a fixed). It follows, in Sec. VE, that 
if we accept current (rough) astrophysical estimates of 
the masses and numbers of inspiralling stellar-mass black 
holes and if we employ only the current LISA design, we 
expect LISA will not see any transitions from inspiral to 
plunge during its lifetime, though it may come close. 

Slight changes in LISA could make some transitions de¬ 


tectable. Dramatic improvements would be required to 
render LISA sensitive to prograde inspirals of stellar-mass 
black holes into rapidly-spinning (a > 0.9) supermassive 
holes. But assuming such inspirals are a small propor¬ 
tion of all inspirals, if the LISA noise curve is lowered 
by a factor 3 (as is under currently discussion for other 
reasons), or if nature provides black holes more massive 
than 10Mg (say 30 Mq) in numbers approaching current 
estimates for 10M e , LISA would have a good chance of 
seeing one or two transitions sometime during its lifetime. 


II. PHYSICAL FRAMEWORK UNDERLYING 
THE TRANSITION LENGTH ESTIMATE 

In the (formal) absence of radiation reaction, a par¬ 
ticle in equatorial orbit about a Kerr hole moves along 
a geodesic. Its radial motion can be determined from a 
first integral of the geodesic equation (equivalent to con¬ 
servation of rest mass; see comments in Appendix [^J: 


dr\ 2 

-J +V[r(r),E,L] = 0, 


( 1 ) 


V = -(E 2 - 1) 2 | (A'-a 2 ^ 2 - 1 )) 2 (L-aEf 

( 2 ) 

Here and throughout this paper all all quantities are, 
for simplicity, made dimensionless using the particle’s 
mass to and the hole’s mass M: E =(orbital energy)/m, 
L =(orbital angular momentum)/mM, r =(orbital 
boyer-lindquist radius)/M, r = (particle’s proper 
time)/M, and a = (hole spin angular momentum)/M 2 . 
Physical solutions may be specified by ( E,L ) or by any 
other pair of equivalent orbital parameters. It is con¬ 
ventional in the inspiral literature to employ as alter¬ 
natives the parameters p (a relativistic generalization of 
semi-latus rectum) and e (a relativistic generalization of 
orbital eccentricity) |Ti], [ll); theseparameters are dis¬ 
cussed in more detail in Appendix [B . 

We concern ourselves with a region of parameter space 
for which the maximum U max of the potential is nearly 
0 (Fig. Q) and which therefore nearly admits a circular 
orbit at the radius r max of the maximum. The geodesic 
equation Eq. ( 0 ) implies that particles can spend an ex¬ 
tremely (logarithmically) long time near the maximum; 
i.e. the particle can “whirl” several times about the hole 
in angle without moving significantly in r. It is conven¬ 
tional to call this portion of the orbit the “whirl.” 

In the presence of radiation reaction (henceforth as¬ 
sumed weak), we must add to the geodesic equation 
(gauge-dependent) time-varying terms which reflect the 
(gauge-dependent) influence of gravitational radiation on 
the test particle’s path. These gauge-dependent terms 
oscillate on the same characteristic timescale as the ra¬ 
diation field. Since the radiation field is predominantly 
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FIG. 1: The effective potential V(r,E,L) for radial geodesic 
motion gradually evolves during inspiral. If the initial eccen¬ 
tricity is nonzero, the maximum of the potential will gradu¬ 
ally decrease until it passes below zero, thus permitting the 
particle to fall into the hole. By way of illustration, we plot 
V(r,E,L) for a = 0 and (solid) (E,L)=( 0.948157, 3.53038) 
and (dashed) ( E,L ) =(0.947454, 3.52092). Under the action 
of radiation reaction, if g = 2 x 10 -3 (an exaggerated mass 
ratio) the first system will evolve into the second after one 
radial period (of the first potential). 


produced during the whirl part of the orbit, the radi¬ 
ation field predominantly oscillates at harmonics of the 
angular frequency Q of circular orbits at the maximum. 
By averaging these reaction forces over a few cycles (e.g., 
over times ~ 2i r/£7) to obtain their secular effect, we 
in principle find expressions for E(t) and L(t). Since 
the averaging time can still (particularly when the parti¬ 
cle whirls several times about the hole) be shorter than 
the time the particle spends whirling around the central 
hole, we can to a good approximation employ Eq. ( 0 ) 
with time-varying E(r), L(t) to follow the orbit when 
the particle is near the maximum of the potential - and 
in particular the interval in which the particle goes from 
nearly-geodesic bound orbit into rapid plunge into the 
hole. 

In this paper, we do not compute E(r) and L{r) i n the 
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thorough, general manner described here. See Sec. 
for a discussion of what information about E(t) and L(t) 
is required for our estimate and how that information is 
obtained. 


A. Why we may still approximate the potential as 
static when computing the radial orbit 

According to Eq. ©> a test particle on an approach to 
a black hole will fall into the hole if the local maximum 
of V, 


I = —V(r max ) = -V m 


( 3 ) 


is negative. Radiation reaction reduces the local maxi¬ 
mum faster than this potential V flattens out. As the 
maximum decreases, the particle spends ever-more of its 
radial cycle near the local maximum. Eventually, we 
reach configurations such as those shown in Fig. [I], where 


the particle can slip over the maximum and fall into the 
black hole. 

While configurations with I ss 0 appear delicately bal¬ 
anced and therefore highly sensitive to small changes in 
{E, L ), in fact under weak conditions (conditions made 
more explicit in Sec. [I I Q ) one may ignore radiation re¬ 
action when computing a radial orbit and treat the po¬ 
tential V as static. 

Suppose a particle starts its whirl with some values 
for ( E,L ) (and therefore I). During the whirl, even 
though (E, L) change, the peak of the potential (/) will 
not change significantly (see Appendix |X]) . Moreover, the 
location of the peak usually moves slowly relative to the 
particle. A co ndit ion for when the latter holds is pre¬ 
sented in Sec. IIC. Therefore, during the whirl one can 
ignore radiation reaction. 

After the particle finishes whirling about the hole, it 
moves outward to its outer radial turning point and back. 
During this period, the maximum does change, from I to 
I' = I + A/. Unless the potential is nearly flat, however, 
the potential away from the neighborhood of the hole will 
not change much as (E, L) change. Agai n, a c ondition for 


when the latter holds is presented in Sec. IIC. Therefore, 


in this interval one can again ignore radiation reaction. 

When we attempt to evolve the particle through the 
next “whirl”, we need the correct value of the height of 
the maximum (now I' = I + A I) to determine how long 
the particle whirls around the hole. Therefore, when we 
start the cycle anew, we must use a potential with pa¬ 
rameters (E+AE,L + AL), with A E and A L the change 
in these constants over the preceeding full radial period. 
If /' < 0, the particle will “bounce” off the maximum 
and we repeat the cycle above once more. But even¬ 
tually we will have I' > 0, at which point the particle 
will move across the maximum during its whirl and will 
subsequently “plunge” into the hole. 

To summarize: so lo ng a s the potential is approxi¬ 
mately static (cf. Sec. [II C| ), we expect we can under¬ 
stand transitions from inspiral to plunge by way of exam¬ 
ining the geodesic equation [Eq. (Q)] in the neighborhood 
of the local maximum, using I £ [0, A I], 


B. Adiabatic approach to estimating the duration 
of the transition from inspiral to plunge 

So long as we can treat the potential as static, we can 
approximate Eq. § in the neighborhood of the poten¬ 
tial’s maximum at r = r max by the form 


7 


d5r 

dt 


Sr 2 


H - y ~ I- 


( 4 ) 


Here 5r = r — r max ; r max is the instantaneously static 
location of the local maximum of V, and also the point 
about which we have expanded the potential; 

r 0 = (y"/2) -1 / 2 


( 5 ) 
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is a constant related to the curvature of the potential at 
the transition location; t is the (dimensionless) time at 
infinity; and 7 is the redshift factor relating proper time 
r to Boyer-linclquist coordinate time at r = r max : 

7 = dt/dT r =r mftx • ( 6 ) 

Estimating the duration of a given transition: Solu¬ 
tions to Eq. ([| give hyperbolic motion; for example, the 
solution appropriate to I > 0 is 

r(t) - r max « V1t q sinht/(r 0 7 ) . (7) 


• Converting to number of cycles: As the parti¬ 
cle passes through the transition region, the par¬ 
ticle “whirls” about the black hole a few times. 
Since its radial location is largely fixed while it 
whirls around the hole, so is its angular frequency 
d(j)/dt = S 2 ; therefore, we can re-express any dura¬ 
tion T c in terms of a “number of orbital cycles” the 
particle “whirls” around the hole N c , defined by 



( 12 ) 


Using this solution, we conclude that the transition time 
going from r - r max = ~Sr lei to r - r max = Sr lei is 


Since we concern ourselves with only Kerr equato¬ 
rial orbits, we have 


T c (Sr re {, I) « 2777 sinh 1 <5r re f/(r 0 v7) 
~ 2777 In 2i5r re f/(r 0 V7) 


( 8 ) 


Hence given (5r re f, a quantity which defines what we mean 
by “the transition extent”, we can estimate the length 
of any transition (characterized by I) at any transition 
location (characterized by the explicit values that go into 
7 , t 0 ). 

Estimating the distribution of transition durations: 
There is no unique transition duration. Rather, we have a 
distribution of durations, depending on the distribution 
of I at the start of the particle’s final whirl. But that 
distribution is simple: since an initial configuration of 
particles will have some distribution of /, since this dis¬ 
tribution evolves smoothly with no “knowledge” of the 
preferred scale A/, and since A I will be smaller than 
any scale in the distribution function, a test particle on 
its final, plunge-triggering whirl has an approximately 
equal probability to have any I £ [0, AJ]. Therefore, the 
probability density for a test particle to have a given du¬ 
ration between T c and T c + dT c is dP oc dT c (dI/dT c ) oc 
dT c exp[—T c / 7 t q ]; see Eq. (||). Denoting by 

Tc- = T c (6r ret , AJ) (9) 


the minimum possible transition duration, and ignoring 
the t iny r egime of transitions which are nonadiabatic (see 
Sec. IIC below), we conclude that 


QM = ijfA. 

r 6 ' 1 + a 


(13) 


• Characteristic duration and variation of T c 
with e: By examining the quadratic approxima¬ 
tion to the potential [Eq. (Q)], or equally well from 
Eq. (J 8 J) , we see that the transition duration is al¬ 
ways T c ~ (few) x 777 — that is, the crossing time is 
around the natural timescale of the effective poten¬ 
tial. Admittedly, since O(AI) = 0(rf), the quantity 
labeled (few) could be — and will be — significant; 
therefore, the logarithmic correction in Eq. (^|) is 
necessary. But for purposes of understanding the 
variation of crossing time with orbital parameters, 
largely we can regard T c ~ 777. For example, we 
expect T c to increase monotonically with decreas¬ 
ing orbital eccentricity e — that is, as the maximum 
possible energy barrier decreases and the potential 
flattens out — simply because 77 does. [By way of 
example, see Eq. ([k|), an expression for 777 appro¬ 
priate to Schwarzchild.] 


• On variation of T c with 77 : Similarly, we can 
loosely characterize the dependence of the duration 
distribution — or, for clarity, T c _ — on 77 by noting 
i) sinh - (2;) « ln 2 x when x is large and ii) A/ oc 77, 
so we can characterize variation with 77 by H(rj 0 ), 
defined by 


dP « 0(T C - T c _ ) e ~( Tc _T <=-)/'t T °dT c / 7 r 0 . (10) 

[where Q(x) =1 when x > 0 , 0 otherwise]. 

We can also characterize distribution of crossing times 
by a function T c {p) such that only a fracti on p of particles 
could (assuming the conditions of Sec. |lIC| hold) have 
longer crossing times. For example, only a fraction 10 -n 
of particles will have duration longer than 

T c , n = T c (Sr leh A I 10" n ) * 2777 In [- ■ 

|_T7l0-"/ 2 \/A7j 

( 11 ) 


Tc-{y) _ x _ In VVo/v 

T c- (VO ) l n ^T-ref/77 -yAfoAJ/T?) 

= ln(?7 0 /?7)JJ(77 0 ) (14) 

[where we have used the fact that A //77 is inde¬ 
pendent of 77 to justify writing the denominator as 
2/ H ( 770 )]. In other words, while the minimum tran¬ 
sition duration will grow slightly shorter with larger 
mass ratios, the dependence (like the dependence 
on AJ/ 77 ) is weak; typically (e.g., for Schwarzchild) 
we find H [0.1,0.4]. 


Additional comments: 















5 


C. Explicit conditions under which we may 
continue to approximate the potential as static 


Throughout our analysis, we have ap proxi mated the 
potential as static. As outlined in Sec. 1IA, there are 


two ways in which this approximation could fail. 

First, the potential away from the maximum could 
change significantly during one whole radial orbit. Gen¬ 
erally the change of V at any specific location is small. 
Such changes therefore matter only if the potential is 
delicately balanced near zero at every point in which the 
particle orbits. More explicitly, we expect problems if the 
change A7 of the potential’s maximum during one whole 
radial orbit is comparable to the difference between the 
maximum and minimum of V. Therefore, we conserva¬ 
tively require 


In lax = v (r max ) - V(r m in) » A/ . 


(15) 


An explicit form for / max is presented in Eq. (BE). Since 
the potential gets very flat as e —> 0, our approximations 
will break down at eccentricities below e m i n , defined by 
solutions to 


A7 — Imaxt^min) • 


(16) 



FIG. 2: The (log of the) change of the potential’s maxi¬ 
mum during the last radial orbit (log 10 [A7/r;]) versus a. The 
points show our AJ for various cases with parameters (p, e) 
close to those consistent with circular orbits (the boundary 
between stable and unstable orbits); these points were ob¬ 
tained from numerical solutions of the Teukolsky equation by 
Gla mpedak is and Kennefic, using the procedure described in 
Sec. IID 1. For each a, solid circles show values of log 10 [A// rj\ 
for several different values of e; that these points are all con¬ 
sistent with a single fit demonstrates that A/is approximately 
independent of e. The solid curve is the quadratic fit Eq. (|l9|). 
[In cases where insufficient data was present for extrapolation 
to the line p = p 3 (e, a), solid circles also indicate rough upper 
and lower limits expected of log 10 [A7/p].] 


Second, the radial location r max of the maximum could 
move significantly while the particle is in its last whirl 
about the hole. Based on Eq. (Q), to prevent against 
this we require I = (dr/dr) 2 !§> (dr max /dr) 2 , i.e. that 

/»w.(^) , = (7*=y. (it) 


The precise procedure that we wi ll us e to estimate 
dr max /dr will be discussed in Sec. IIP . In summary, 
so long as I / a d,min, or equivalently so long as the 
crossing duration significantly shorter than 


Fc, a d,min 


Tc(dr Te f, 7 a d ,min) 


(18) 


gradual motion of the potential will not significantly alter 
the transition length estimates presented earlier. 


D. Inputs necessary for estimating the transition 
length 

In the above we have outlined a computational proce¬ 
dure which takes as input Sr le f and knowledge about ra¬ 
diation reaction (namely, about A7 and about dr max /dt) 
and which gives us in return an estimate of the length 
of any specific transition from inspiral to plunge. We 
now describe the explicit approximations we shall use to 
estimate A7 and 7 a d,min from known information about 
E(t) and L(t). We also make an explicit choice for <5r re f. 


1. Estimating A7 


As described in Sec. IIA, we obtain A I by comparing 
the potential V when the conserved constants are (E, L) 
to the potential V when they are (E + A E,L + AT), 
where A E and A L are the change in the appropriate 
conserved constants over one radial orbit. We obtain 
AE and AT from numerical solutions to the Teukolsky 
equation. From their code, Glampedakis and Kennefick 
have kindly provided time-averaged fluxes ( dE/dt ,) and 
( dL/dt) ]lT| , which, when combined with an expression 
for the radial period T(E,L) as given in any classical 
relativity text pe.g., Eq. (33.37) of MTW ||], yields A E 
and AT, and thus A7 [Eqs. (|j) and (||)]. In this fashion, 
for each black hole (parametrized by spin parameter a), 
we can find AI(p,e,a) for any equitorial geodesic with 
parameters (p,e). 

We need A7 only for the particle’s last whirl. Since the 
maximum is extremely close to zero, the orbital param¬ 
eters ( p, e, a) nearly satisfy a condition for the existence 
of (unstable) circular orbits p = p s (e,a) (see Appendix 
|b|). This curve also necessarily serves as the boundary 
between stable orbits and plunge. In the vicinity of this 
boundary line, A7 is well-approximated by its nonzero 
values on the boundary. So for our computation we seek 
an expression A7(e,a). 

In practice, from the values of A7 at points near this 
boundary line, we extrapolate to estimate A7 on the 
boundary surface itself. Figure || shows the results of 
our extrapolation. 

One can argue that AI(p, e,a) on the last-stable-orbit 
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boundary p = p s (e, a) should largely be independent of 
e at moderate eccentricity Jl^|. For this reason, Figure 
shows results only as a function of one parameter (a). 
Numerical data over the range e £ [0,0.5] support this 
conjecture. Therefore, so long as we avoid e « 1, where 
this conjecture has not yet been tested and likely fails, 
we can approximate A //77 by a function independent of 
e. Fitting a relatively simple function (exponential form 
in a, independent of e) to the data in Fig. we find 

log 10 — = -0.8972 + 0.7911a + 0.3047a 2 . (19) 

V 


2. Estimating I a d,n 


To evaluate I a d,min, we need no more than i) knowledge 
of the potential (which tells us r max as a function of E, 
L) and ii) knowledge of dE/dt , dL/dt when the particles 
are in nearly-circular orbit near the hole. 

In principle, we could approximate the latter by the 
appropriate values for an exactly circular (unstable) or¬ 
bit. As a practical matter, comprehensive tabulation 
of the physically appropriate instantaneous dE/dt and 
dL/dt for all transitions of interest — namely, the values 
appropriate to a circular unstable orbit — proves time- 
consuming and technically challenging. Furthermore, be¬ 
cause the crossing time depends only weakly (logarith¬ 
mically) on /, and because exceedingly few particles will 
have I < 7 a d,min 5 we only need / a d,min to order of magni¬ 
tude. 

Therefore, for practical purposes, when estimating 
Aid,min by way of Eq. (|l7]) we will i) perform the compu¬ 
tation for (ir ma x/dT analytically in terms of dE/dt and 
dL/dt , ii) simplify under the assumption dE = LldL, 
which would be valid if we used the true forms for 
dE/dt and dL/dt, and then iii) insert for dL/dt the 
Peters-Mathews expression (an estimate obtained using 
linearized, quadrupolar emission from newtonian orbits) 


with circular orbits there. If we were to try to detect 
these gravitational waves — for simplicity, focusing on 
the dominant frequency component, to = 2 f 1 — we would 
want to insure that our model </> = 2Qt for the gravita¬ 
tional wave phase agrees, within n, with the true wave 
phase. 

The true rate of change of orbital phase is 


m r \ 

H (r ’ E ' L) 


g^L—g^E 
—g u E + gt't’L ’ 


( 21 ) 


where g a P are known Kerr metric functions in Boyer- 
Lindquist coordinates, and E, L are consistent with the 
circular orbit at r = r max (use standard expressions for 
E, L appropriate to circular orbits, such as Eqs. (2.12) 
and (2.13) of Bardeen, Press, and Teukolsky Q). De¬ 
manding that the difference between the true angular 
phase and our fiducial reference 2 f It be no more than 7 r 
over the length of the transition, we find a constraint on 
the crossing duration T c : 



— (r )' 
dt [ max) 


dt. 


( 22 ) 


When we insert r(t ) — r max = Asinh[t/r 0 7 ] with A = 
VIr 0 [Eq. (Q)] into the above, we find an expression we 
can invert for Sr re {(AI ): 


7t/4 


d(d(j)/dt) 

dr 



d{d(j)/dt) 

dr 


jToSr re f 


(23) 

[where the constant A has been neglected in this expres¬ 
sion, as it is always much smaller than <5r re f]. Solving for 
Sr le f, we obtain 


Ar r ef 


7T 

47t 0 I d{d<j)/dt)/dr \ 


(24) 


dL 

dt 


32 1 


(1 - e 2 ) 3 / 2 


l + le 



We will use this form even when it predicts <5r re f = 0(1) 
[in other words, when 8r re f = 0(M ) when we convert to 
physical units]. Notice this Sr le f is independent of mass 
ratio. 


3. Choosing 5r re f 

To complete our procedure, we must define “the” tran¬ 
sition duration. Unfortunately, because “the” transition 
from inspiral to plunge occurs at no definite location, has 
no well-defined start or finish, the transition duration re¬ 
mains a matter of convention [jn|. We shall adopt a 
convention motivated by a simple model of gravitational- 
wave data analysis. 

The key feature of waves emitted during the transi¬ 
tion is their considerable simplicity: they are emitted 
from a nearly-circular-equatorial orbit at r max , and hence 
are characterized by the angular frequency U associated 


III. ESTIMATING THE PROBABILITY FOR 
LISA TO OBSERVE A TRANSITION 

We wish to estimate, for each choice of the supermas- 
sive hole’s angular momentum and distance from earth, 
and for each choice of test particle orbital parameters, 
the signal-to-noise ( S/N ) LISA would obtain from waves 
emitted during the transition. By combining this S/N 
with the (poorly-known) statistics of black-hole inspirals, 
we can estimate the probability LISA will see a transition 
event (e.g., have S/N > 5). 
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A. Estimating LISA’s signal to noise for a given 
transition 


Since the transition waves are emitted by a circular 
orbit of frequency 

/orb = f2( T* max )/2 t tM, (25) 

the gravitational waves will be at that frequency and its 
harmonics. For simplicity, assume that LISA detects only 
the strongest waves, the waves emitted from the second 
harmonic /t r = 2/ or b- These waves will last for an inter¬ 
val 

At = N c // orb (= MT c ) . (26) 

We can approximate their characteristic rms (source- 
orientation-averaged) amplitude [following OT equation 
(4.7)] as a Peters-Mathews-style quadrapole term (aver¬ 
aged over all orientations) times a relativistic correction: 


h lms 


8 M 
7 t~d 


7yl7(rmax) 


2/3 



(27) 


Here d is the distance to the source and £ 00,2 is a rela¬ 
tivistic correction factor defined explicitly in OT equation 
(2.3). 

LISA has a spectral density of noise Sh for waves inci¬ 
dent on it with optimal propagation direction and polar¬ 
ization; it has spectral density 5 Sh for typical directions 
and polarizations. Therefore, on average, LISA should 
accumulate a signal-to-noise from the transition event 
given by 


(S/N) tr 


/i rms 

V5S h (/ tr )/Ai ' 


(28) 


Particularly special sources could have significantly 
higher S/N. For example, we can pick up an increase 
of v/5 if the source is ideally positioned on the sky, and a 
similar increase if the source itself is optimally oriented. 
But overall, the above scheme suffices to estimate the 
signal-to-noise LISA would see from the transition be¬ 
tween inspiral and plunge for any capture m into M with 
any specific source parameters (e.g., e, a) at any distance 
d. 


Explicit expressions needed to compute LISA’s signal to 
noise for a given transition 


To evaluate Eq. (|2§|), we need in addition to N c and 
H [which enter into S/N via St and / tr ] the LISA noise 
curve Sh and the relativistic correction factor £ 00 , 2 - The 
LISA noise curve may be modeled by [OT equation (4.9)] 


Sh(f) 


(4.6 x 1CT 21 ) 2 + (3.5 x 10~ 26 ) 2 (yy^j 

+ < 3 - 5 * 1(r,, > 2 (t k 



The appropriate relativistic correction factor £ 00,2 can in 
principle be extracted from simulations of waves emitted 
by particles in unstable circular orbits. As in practice the 
latter proves time-consuming to evaluate and tabulate for 
all possible eccentric orbits and for all a, for simplicity 
we will assume that the appropriate relativistic correc¬ 
tion factor is i) fixed for all orbits close to a black hole 
of angular momentum a and ii) given explicitly by the 
value appropriate to the innermost stable circular orbit 
(ISCO). This latter expression has been tabulated by Ori 
and Thorne (see the £ 00,2 column in their Table II); we 
approximate their results by 

log 10 (£ 00 , 2 ) qt « -0.0473 + 0.211a: - 0.053a; 2 

+0.034a; 3 + 0.010a: 4 (30) 

where x = log 10 (l — a). 


Dominant terms in the signal-to-noise estimate 


As written, the signal-to-noise estimate Eq. ( p8| ) dis¬ 
guises what kinds of effects predominantly influence it 
for example, whether changes in the strength of radi¬ 
ation emitted prove more important or less than changes 
in the duration T c of the transition. To clarify the domi¬ 
nant contributions to our estimate, fix some a and com¬ 
pare the signal-to-noise between two transitions (1,2) in¬ 
volving otherwise arbitrary parameters (e.g., ?n, M, d , 
e). Substituting expressions for At [Eq. (|2^) [, /i rms [Eq. 
(§ 3 )], and / tr = 2/ orb [Eq. (§§)] into Eq. (@); assuming 
£ 00,2 is a fixed function of a; and comparing the resulting 
S/N at two sets of orbital parameters, we find 


(S/N) 2 = m2di [N^ / Hi \ 1/6 

(s/n), m,d 2 y iv Ci i v^y 


M,S h 

( nM 2 ) 

M 2 S h 



(31) 


The first two terms reflect the natural m/d scaling of 
emitted waves. The third term reflects the fact that more 
orbits around the hole during the transition mean more 
gravitational wave cycles seen by LISA. The fourth term, 
which combines the fact that gravitational waves emitted 
closer to the hole are stronger and yet last for less time, 
is to a good approximation constant. Finally, the last 
term reflects LISA’s sensitivity. The only term which 
depends explicitly on M (ignoring the weak variation in 
N c ), this last term selects black hole masses which have 
their transition close to optimally positioned in the LISA 
band, or M ss (few) x 10 6 . So long as the mass is so, this 
term varies comparatively little. 













B. Method for estimating the probability of 
detecting some transition during LISA’s operation 


Above we gave a procedure for computing the S/N 
for any given source. But the sources which produce the 
strongest signals (inspirals very close by) are rare. There¬ 
fore, for any given ( S/N) 0 we have a certain probability 
that, during the entire operation time Tj J of LISA, we 
detect no inspirals with S/N > ( S/N) a . 

Since the relevant statistics for supermassive black 
holes and compact objects are poorly known, we will not 
attempt a detailed calculation that allows for all possi¬ 
ble factors (e.g., source-orientation effects). Instead, for 
a first-pass estimate of the likelihood that LISA will see 
a transition, we will i) fix M = 10 6 , ii) approximate 
LISA’s noise curve as flat (in other words, ignore varia¬ 
tions in S/N due to the emitted radiation being slightly 
off LISA’s peak sensitivity), iii) ignore any orientation- 
related increase in the emissivity of the source or the 
sensitivity of LISA, iv) approximate N c as independent 
of to, v) further replace N c at each a by some characteris¬ 
tic number of cycles (the precise value to be chosen later, 
when we understand how N c varies), and vi) assume all 
black holes have the same value of a (again, to be chosen 
later). To be particularly explicit, we assume the S/N 
varies with to, d, and a in the following manner: 



(to, d, a) 


to lGpc / S \ 
1OM 0 d \n) a 


(32) 


Here (S/N)a = (S/N)a{ 10M q , lGpc, a) is a fiducial ap¬ 
proximation to the signal to noise ratio for an inspiral 
with to = IOMq, d = lGpc, and a. 

Suppose we have a discrete family of possible com¬ 
pact objects of masses m.k with rates (per galaxy con¬ 
taining a 1O 6 M 0 hole) suppose the number density 
of galaxies containing a 1O 6 M 0 hole is p g . Subdivid¬ 
ing the universe into cubes of cell size Ar, we find the 
probability a given cell has an inspiral of mass nik into a 
1O 6 M 0 hole at some time during the lifetime T& of LISA 
is Pk = PgfkTLNr 3, . Suppose we’re concerned with a 
threshold S/N level S/N = s a . At such a level we could 
see a source of mass to*, out to a distance dk = Krrik/s a . 
If no inspirals have S/N > s 0 , then for every cell in 
range, we have no inspirals of any mass type. Therefore, 
the probability that no inspirals occur with S/N < s 0 is 


P(no S/N > s 0 ) 


nu-*) 4 "^ 


exp 


AnR net T L K 3 (to 3 ) 


(33) 


where in the last line we use pk 1, R ne t = Y2,Pg r k 
(the net event rate per unit volume for all inspirals), and 
< to 3 >= p g rkm 3 /R ne t (the mean cubed mass of in- 
spiralling bodies, where weights are by event rate). Note 


that 4ttI\ 3 ( to 3 ) /(3s 3 ) is the volume of space in which 

1 /3 

an inspiral involving a mass (to 3 ) can be seen with 
a signal-to-noise ratio > s D . Necessarily, the probability 
that some source has S/N > s 0 is P(some S/N > s 0 ) = 

1 — P(no S/N > s 0 ) 

We can reorganize this expression to tell us, for a given 
probability P no , at what S/N we will have a probability 
P no of having no signals of stronger strength: 

(S\ (P s = (m 3 ) 1/3 lGpc _1 

W/J noJ ■ 10M 0 { 4«R netTL) -V 3 Ml/P no )]i/ 3 

-(«, 


C. Probability of detecting a transition during 
LISA’s operation 


The S/N threshold [Eq. (34)] depends very sensitively 

1 /3 

(through (to 3 ) ) on low-probability high-mass inspi¬ 

rals. By way of illustration, a family of O.6M 0 white 
dwarfs inspiralling with rate R and a black hole family 
of mass 30M 0 and rate 10 _4 P contribute in similar pro- 

portions to (to 3 ) . At present, the astrophysical com¬ 

munity lacks sufficiently understanding of the high-mass 

tail to be able to reliably compute (to 3 ) . Therefore, 

we will neglect such objects and focus on the slightly- 
better understood problem of capture of conventional 
compact objects. Doing so, we will underestimate the 
true ( S/N) no . 

Even disregarding the high-mass holes, event rates for 
capture [O, p~5| remain very loosely determined, ranging 
from rates of ~ 2 x 10 _6 /yr/galaxy to ~ 10" 4 /yr/galaxy. 
We take two cases as characteristic: 


• Freitag (F) Based on astrophysical discussion by 
Miralda-Escude and Gould [[Ig], Freitag allows for 
three species: white dwarfs (towd = 0.6M 0 , 
rwD ~ 10 _5 /yr); neutron stars (tons = 1-4M 0 , 
r/vs ~ 2 x 10~ 6 /yr); and black holes (tobh = 7AP 0 , 
tbh ~ 10” 6 /yr). In this case, the net event rate 
P ne t is dominated by low-mass WD inspirals, but 
black holes dominate the events seen by LISA. 
Using a LISA lifetime Tj, = 3 yr and (based on 
Sigurdsson and Rees’s estimate that the density 
of 1O 6 M 0 holes at their cores is around the den¬ 
sity of spirals, since spirals have low mass and el¬ 
lipticals high mass supermassive holes jn|) p g ~ 
0.003/Mpc 3 , we find 


N 


9 «7 ( ln ( 2 ) N V3 

■ 67 U(l/P„o), 


N 


(1OM 0 , lGpc). 
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• Sigurdsson and Rees (SR) They consider two types 
of galaxies — spirals and dwarf ellipticals — but 
only the latter leads to significant event rates. In 
that case, they uses the following masses and rates 
for the three species: WD (towd = 0.6M 0 , rwD ~ 
3 x 10 _8 /yr); NS (tons = 1.4M 0 , tns ~ 10~ 7 /yr); 
and BH (tobh = 5M 0 , tbh ~ 10 _6 /yr). (For 
black holes, these authors provide only an off-the- 
cuff estimate; we have taken some liberty in inter¬ 
preting it, choosing a mildly optimistic character¬ 
istic black hole mass.) Again using = 3 yr and 
p g = 0.003/Mpc 3 , we find 

m =igo ( M2) y /a 

U;„o,SR ' Vln(l/Pno 

x (1OM 0 , lGpc). 


In performing both calculations, we use Sigurdsson and 
Rees’s estimate that the density of galaxies with (spirals 
have low-mass holes; ellipticals and others tend to have 
more), so p g ~ 0.003/Mpc 3 . Also, we use a LISA lifetime 
T l = 3 yr. 

In sum, we suspect that on astrophysical grounds we 
will have a 50% chance of seeing no source with S/N 
roughly greater than 


2.5 x 



(35) 


where the (S/N )a will chosen to be the most reasonable 
S/N over all orbital parameters (e) and black hole spins 
(a), given the fiducial parameters d = lGpc, m = 1OM 0 , 
and M = 10 6 M 0 . 


IV. SCHWARZCHILD SUPERMASSIVE BLACK 
HOLE (SMBH) 

To illustrate this scheme in a case where all terms are 
algebraically tractable, we discuss the range of probable 
transition durations when the capturing hole has no an¬ 
gular momentum (Schwarzchild). 


Since we have only two free parameters, the three roots 
are not independent; they satisfy a self-consistency poly¬ 
nomial. For this reason, we introduce p,e — parameters 
analogous to semi-latus rectum and eccentricity in clas¬ 
sical mecha nics . Employing a consistency relation [gen¬ 
erally Eq. (B 6 ) of Appendix (b| to set f , we find 


r± = 


P 


1 T e 


J2p_ 
p — 4’ 


(37) 


F 2 = ip - 2 ~ 2e) (p - 2 + 2e) 2 = p 2 M 2 

p(p- 3-e 2 ) ’ p- 3 —e 2 ' 

(38) 

These p, e parameters have the notable advantage that 
bound orbits (orbits that cannot escape to infinity) and 
non-plunging orbits (orbits which avoid the central black 
hole) are easy to describe: bound orbits have e £ [ 0 , 1 ], 
while non-plunging orbits have 0 < r_ — f = p(p — 6 — 
2e)/[(l + e)(p — 4)], or 

z = p — 6 — 2e>0. (39) 

As one approaches the transition, the maximum of the 
potential V decreases, r_ approaches f, and z approaches 
0. We can equivalently specify the location r of a tran¬ 
sition by only one of p or e, with the other determined 
by z = 0. I will typically use e. For example, a tran¬ 
sition of eccentricity e occurs at radius r = r_ = f = 
2(3 + e)/(l + e). 

The parameters p, e used here are identical to those 
used the Teukolsky-equation-based inspiral literature 
® 0 - For example, the above discussion mirrors 
that in Cutler, Kennefick, and Poisson 0 between their 
equations (2.4) and (2.8), with the change of notation 
r\ —► r_, T 2 —> r+ and r 3 —+ f. 


B. Dependence of Transition Parameters on 
Eccentricity 


We know the potential [Eq. (36)]; hence we find that 
when r_ = r (=the transition radius) we have 


A. Choosing parameters 


V" = -(1 - E 2 )2^—^— (40) 


Rather than using E , L to characterize the orbit, when 
the orbit is confined in radius between two turning points 
(i.e. when it is bound), it is far simpler to characterize 
the potential V = —(dr/dr) 2 by the location of its 3 
roots, r ±, f , where r± are the inner and outer turning 
points of the bound orbit and r is the innermost root: 

■2 

— (r_|_ — r) (r — r_) (r — r) 

= e 2 -(H)( i+ £) < 36 > 



and therefore, substituting r± and E from Eqs. (0), 
( 0 ) into t 0 = (V”/2)~ 1/2 we find 


To 


(3 + e) 


2(9 — e 2 ) 
e(l + e ) 3 


(41) 


Similarly, substituting r = r_ into 7 = — g u E gives us 
the characteristic time required to make the transition: 


2(3+ e ) 2 
y 7 e(l + e ) 3 


(42) 
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Since d(j)/dt 
find 


g^L/{-g tt E) 


L/r 2 (l — 2/r)E, we 


f d_ df\ 

\dr dt) 


T —^max 


2 (r 


5/2 , 

T max y max 


3)_ 

- 2 ) 


(43) 


when E, L are consistent with a circular orbit at r = 
r m ax- Using the above and Eq. (|24|), we conclude that 
for moderate eccentricity the natural “transition extent” 
<5r ref is 


3c re f 


\/2e(3 + e) 
7r (l + e)(3-e)' 


(44) 


This scale naturally varies with the scale of the potential 
(namely, Sr le f oc 1/y/V") as e —> 0. 

Further straightforward computations using the poten¬ 
tial [Eq. ([36])] reveal how r max varies due to loss of E and 
L via wave emission when the particle is nearly in a cir¬ 
cular orbit (so dE « QdL): 

dv max dL ^ (r Tn ax 3) 

dt dt r max - 6 


We therefore can express / a d,min in terms of (tab¬ 
ulated) known radiation-reaction angular momentum 
fluxes ( dL/dt): 


T t 'dL\ 2 (3-e) 2 (3 + e) 

ad ’ min "" 1 dt J 2 e 2 (l + e) 


(45) 


To obtain a rough approximation of / a d,min, rather than 
use the true dL/dt appropriate to circular orbits, we ap¬ 
proximate dL/dt by the Peters-Mathews expression [Eq. 

©]■ 

Our scheme ceases to apply when the eccentricity is 
below e ad defined by A I = / max (e a d) [Eq. (@]. In the 
special case of a = 0, where p s = 6 + 2e, the definition of 


I max [Eq. ( |B8|) 1 reduces to 


Imax — 


32 


27 (9 — e 2 )(l + e) 


(46) 


C. Transition duration 



FIG. 3: Plots of various predictions for the expected num¬ 
ber of angular cycles (N c ) versus eccentricity (e) for a tran¬ 
sition from inspiral to plunge with g = ICG 5 . The top solid 
curve is the number of cycles when T c = T/ a a [Eq. (|lq)], an 
estimate of the longest possible (adiabatic) transition dura¬ 
tion. The bottom solid curve is the number of cycles when 
T c = T c [Eq. (||)]. The 6 curves in between are the number 
of cycles when T c = T 0 ,i, . ..T c ,6 [Eq. (Ill)]; as only a frac¬ 
tion Rs 10 -1 ,..., 10 -6 of inspirals can have durations above 
these curves (respectively), they illustrate how few particles 
have durations significantly differing from T Cj _. The dot in¬ 
dicates the Ori-Tliorne (OT) prediction for circular equato¬ 
rial inspiral. The dashed curve is a characteristic-scale-based 
prediction based on T c = 47t 0 , used to illustrate the sig¬ 
nificance of the “logarithmic correction.” The plot starts at 
e = e a d ~ 0.0215, at which point A/ = J max . 


substantially less. At low eccentricity, most transitions 
seem to approach a result somewhat different than the 
OT circular estimate. Since OT use a different conven¬ 
tion for Sr [^0| and since significant changes could still 
occur in the fundamentally nonadiabatic region between 
e = 0 and e = e a d, we do not find the discrepancy trou¬ 
bling. 

In the above, we show results for only g = 10~ 5 (say, 
for m = 10 Mq and M = 10 6 ). As the variation of the 
duration with g is weak — we find = 10~ 5 ) € 

[0.1,0.4]) [Eq. @)] — even substantially different test 
particle masses (e.g., m £ [O.1,3O]M 0 with M = 10 6 ) 
lead to results of the form above, scaled up or down by 
a factor < 2. 


With all the necessary elements assembled, we can ap¬ 
ply our program [Eqs. (g), ©, @] to estimate the 
distribution in number of orbital cycles N c = T c Ll/2n we 
expect when a particle spirals into a nonspinning black 
hole at some fixed, known eccentricity e. 

The results for g = 10~~ 5 are shown in Fig. g. When 
our adiabatic approximation applies, we find that to a 
good approximation (within around 1 cycle) most tran¬ 
sitions should have duration close to the shortest transi¬ 
tion duration _/V c ,ad = T Cja dU/27r £ [3,5]. In particular, 
within the region e > e a d where our adiabatic approxima¬ 
tion applies, almost all transitions will last for less than 
the Ori-Thorne (OT) circular duration; most will last for 


D. Prospects for LISA detecting a given transition 


As discussed in Sec. HI, we can estimate the signal- 

11 ) • 


to-noise ratio for a given transition using Eq. 

For the standard case of a 10Mq particle falling into 
a M = 10 6 Mq, application of that formula reveals no 
higher S/N than that predicted by Ori and Thorne; 
moreover, barring astrophysically unlikely masses, all 
transitions have too low a S/N to be detected. [See 
Fig. 0 below for details.] For example, if an inspiral 
of mass to into a 10 6 Mq hole occurs at the fiducial dis¬ 
tance 1 Gpc with e = 1/3, we have a 90% chance that 
S/N £ [0.91, 1.01](to/ 10M q ). 
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The results for S/N can be well-approximated by way 
of Eq. ©) and a comparison with Ori and Thorne’s 
results for circular inspiral. (See Appendix ^ for a sum¬ 
mary of OT results). Specifically, using the fiducial 
case of 1OM 0 on 1O 6 M 0 at 1 Gpc, for which we have 
( S/N)ot = 1-6 and N c ,ot = 10.5, we find the general 
S/N for captures by a M = 1O 6 M 0 hole to be about 


(S/N) 


1.6 x 



to lGpc 
1OM 0 d 


V. KERR SMBH 

The Kerr case follows similarly, save with an additional 
parameter (a). 


A. Parameterizing Orbits, Potential 


As before, it is simplest to characterize the potential 
by its three roots r± = p/(l T e), f: 

( dr \ 2 l - E 2 

=-V= (r + -r)(r- r_) (r-f). (47) 

And as before we can define r± =p/(lye); as before, 
we find a self-consistency relation P(p, e, a,f) [Eq. (B 6 )], 
permitting us to solve for f(p,e 1 a). As before, we can 
characterize the proximity to the last-stable surface by 
way of the separation between the two innermost roots 
(r_ — f). Finally, as before, for each fixed black hole 
(a =const) and each particle exactly on the transition line 
from orbit to plunge, the particle can have e £ [ 0 , 1 ); its 
p will be constrained by the analogue of the Schwarzchild 
p = 6 + 2e: the self-consistency relation Eq. (B7), which 
implicitly defines p s (e, a ) such that r(p s , e, a) = p s /( 1 + 
e). 


B. Dependence of transition parameters on e,p 


Since the potential has the same structure as before, 
the same general expression Eq. (f!o|) applies, with E now 
determined by expressions in Appendix |Bj. By explicitly 
differentiating the potential [Eq. (||)], inserting the defi¬ 
nitions r± = p/( 1 =F e), and demanding the inner turning 
point is a maximum (so r = r_ = p/( 1 + e)), we find 


V" = - 8 e 


(1 + e ) 3 
(3 - e)p 3 


(48) 


and therefore know r Q = ( V"/2)~ 1 / 2 in terms of p,e at 
the transition. 

The 7 factor follows from inserting r = p/( 1 + e) into 
the usual expression for dt/dr : 


7 =| : : = -g a E + g^L 


(49) 


where g tt and g^ are known Kerr metric coeffici ents . 
Here, E and L are evaluated using the expressions (B4) 
and ( [ESS] ) discussed in Appendix [j|, with f = p s /(l + e). 

We obtain the transition extent <5?' re f with the usual 
Eq. (p4|). This requires 7 [Eq. (f49|), above], r D (also 
above), and d(d<j>/dt)/dr [Eq. ©)]. 

Finally, as in the Schwarzchild case we estimate / a d,min 
[Eq. ([[ 7 ])] and thus T c>a d [Eq. ©)] via i) expressing 
dr max /dr in terms of dL/dt using explicit expressions for 
t max (P , L) and dE = fldL , giving 


dr max dL 2r(2a v / r + r 2 — 3 r ) 3 / 2 
dt dt (r 3 / 2 + a)(r 2 — 6r + 8a^/r ~ 3ci 2 ) r=r ’ 

”(50) 

(where we have used the orbital parameters E,L consis¬ 
tent with a circular orbit at r = r max Jfqj); then ii) using 
the Peters-Mathews expressions for dL/dt [Eq. (g^)] to 
construct an approximate expression for dr max /c?t, which 
we then iii) insert in Eq. ( |f7| ) to estimate the boundary 
between adiabatic and nonadiabatic transitions. 


C. Transition duration 


Combining these together, we can deduce the range 
of plausible transition durations for a test particle of ec¬ 
centricity e falling into a hole of angular momentum a, 
measured as number of orbital cycles N c (e, a) = T c Ll/2ir. 
Plots of the number of cycles appropriate to T c = T c _ 
[Eq. (§)], to T c = T c>ad Eq. (|T|)], and to T c = T c> 1 [Eq. 
(|il|)] appear in Figs. 11 and|^, respectively. 

These plots all assume a fiducial source (fOM 0 on 
fO 6 M 0 ). In these plots, we truncate the range of e, a 
allowed because, i) we need e larger than e a d [Eq. ©]; 
and ii) realistic astrophysical black holes have a < 0.998 
JTt| . Also, in these plots, we do not extend to e « 1 
because we do not have data for I in this region, nor do 
we expect our estimate of A I [Fig. [|] to be reliable in 
this extreme. 

At each a, we see behavior largely similar to the 
Schwarzchild results discussed in Sec. |v|: i) almost all 
transitions take less time than the Ori-Thorne result for 
e = 0 ; ii) as we increase the eccentricity, the transition 
duration decreases; and iii) since T Ci 1 « T c> _ (compare 
Figs. [I] ancl|(i|), most transitions last close to the shortest- 
possible transition duration. 


D. Prospects for LISA detecting a given transition 

As in the Schwarzchild case, since eccentric usually 
transitions last for fewer angular cycles than their circular 
analogues, they are less detectable as well. Thus, in the 
fiducial case of captures of a 1OM 0 hole by a IO 6 M 0 hole, 
the data from Ori-Thorne Table II provides an upper 
bound on the S/N seen by LISA (shown in Fig. jfj). Since 
this bound is small, we have little chance of seeing any 
given transition. 
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FIG. 4: These two plots illustrate the shortest possible num¬ 
ber of cycles N c = £IT c -/2t: [Eq. (^)] a transition could last, 
versus eccentricity (e) and black hole angular momentum (a) 
for the fiducial source (IOMq into 10 6 Mq). In both plots, con¬ 
tours are cut off, and bounding curves appear (shown heavy 
solid), when A7 = 7 a d, m m and when a = 0.998. 

One should notice, however, that the distribution of 
S/N with orbital parameters is very flat and not much 
below 1. Therefore, only a modest improvement in 
LISA’s noise spectrum Sh could render most (measured 
by volume of parameter space) of the transitions de¬ 
tectable. 


E. On probability of detection 

Because LISA at present has so poor prospects for de¬ 
tecting the “fiducial” source (to = 1OM 0 at lGpc), it has 




FIG. 5: These two plots illustrate the longest possible number 
of cycles an (adiabatic) transition could last N c = f2T Cja d/27r 
[Eq. (|l^)1, versus eccentricity (e) and black hole angular 
momentum (a) for the fiducial source (10A7g into 10 6 Mg). 
Transitions of such long duration are extremely unlikely un¬ 
less e « e a d [see Eq. (|l6|), (|f(ij)]. In both plots, contours 
are cut off, and bounding curves appear (shown heavy solid), 
when A7 = 7 a d, m in and when a = 0.998. 


a poor chance of seeing any source at all. Even assuming 
all LISA sources had orbital parameters chosen to give 
the longest-plausible transition length (the OT circular 
inspiral duration, which has S/Na < 1-6), by the esti¬ 
mate of Eq. ( |35| ) we expect we have a ~50% chance of 
no signal with S/N > 4 being present in the datastream. 
With more realistic orbital parameters, we would expect 
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FIG. 6: These two plots illustrate N c = QT c ,i/2n versus e 
and a for the fiducial case. When the results of Fig. [s| are 
appreciably larger than these — that is, everywhere except 
near the left boundary e = e a d — ~ 90% of transitions with 
orbital parameters e, a will take less than the number of cycles 
shown in this plot to traverse the transition region. 


a ~50% chance of no signal > 2.3. In other words, LISA 
has a good to excellent chance of not seeing any transi¬ 
tions from inspiral to plunge in its lifetime. 

A modest improvement in LISA’s noise curve, how¬ 
ever, would make a few circular (and to a lesser degree 
eccentric) transitions from inspiral to plunge detectable. 



FIG. 7: This plot illustrates S/N [Eq. (pj|)1, given a fiducial 
source (10Mg into 10®Mg, at lGpc) given that the black hole 
has angular momentum a, the transition occurs at eccentricity 
e, and given the transition takes the shortest possible time (or 
T c = T c ,~) [Eq. ([)])]. As argued in the text, this time (and 
thus this S/N) will be to a good approximation characteristic 
of all transitions with those interaction parameters. 


VI. SUMMARY 

This paper has introduced a framework (depending on 
observational or other conventions) that extends the Ori- 
Thorne prediction for the transition duration from in¬ 
spiral to plunge to include eccentric orbits. While the 
framework and applications contain many oversimplifi¬ 
cations — most notably, the fit to A I(e,a) and and the 
lack of a physically meaningful convention for Sr Te { — the 
essential physics should be captured by Sec. jn| 

This paper then applies that framework to probable 
LISA sources to suggest that, because an eccentric tran¬ 
sition is generally only slightly briefer than a circular one, 
LISA should have only slightly worse prospects to resolve 
the transition from inspiral to plunge for eccentric orbits 
than for circular ones. While the prospects for detecting 
circular (and hence eccentric) transitions with LISA are 
not good, they are not necessarily bad: modest changes 
to the LISA noise floor could render a signal marginally 
detectable. Therefore, more detailed investigations could 
be of use. 

Potentially, we could use other portions of orbits that 
pass close to the hole for example, the previous few 
“bounces” off the inner portion of the radial potential - 
as probes of the strong-field metric. Analyzed separately 
(using the same A I framework) each of these “bounces” 
should provide in itself at best of order the same S/N as 
the transition. If the source has already been detected 
with good confidence, we should be able to coherently 
integrate over many such bounces and build up excellent 
S/N. 
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Finally, we could hope that eccentric inclined orbits 
might, by some happenstance of parameters, admit a 
regime of significantly longer transition times. The 
prospect seems unlikely, but the author may address it 
in a future paper. 


APPENDIX A: EVOLUTION OF THE MAXIMUM 

As the conserved constants E , L evolve, the height 
of any local maximum in the potential V [Eq. ([![)] will 
similarly evolve. Because we are at a local maximum, we 
can find a simple expression for the rate of change of the 
value of the potential at that local maximum: 

U max (£, L) = V(r max (E, L), E, L ) (Al) 

(where r m ax, the location of the potential’s local maxi¬ 
mum, is a solution to dV/dr = 0). In general, to order 
of magnitude, we expect it goes as ~ U/T gw for r gw the 
gravitational wave timescale r gw ~ E/E ~ L/L. But 
when the particle is nearly on a circular orbit, then the 
source of radiation nearly satisfies helical symmetry and 
therefore dE « f IdL for SI the angular frequency of the 
circular orbit. And in these special conditions V changes 
even more slowly than we would normally expect. 

To be explicit, we evaluate dV max /dt , which (because 
we are at a maximum) we can generally express as fol¬ 
lows: 


dU max dV IdL i dEdV/dE 
dt ~ ~dL \~dt + ~dt dV/dL 


(A2) 


We can most transparently prove the necessary result 
by rewriting the potential V more abstractly than the 
standard form of Eq. (Q). Recall one derives the ra¬ 
dial potential from the constancy of the test-particle’s 
rest mass [e.g., g ab p a Pb = — m 2 }. Since the Kerr metric 
in boyer-lindquist coordinates has form g a g = gttdt 2 + 
2 gt^dtdcj) + g^dc/ 2 + g rr dr 2 + gggdd 2 , by employing the 
definitions of E and L (e.g., —pt = energy) and the 
existence of an equatorial orbit, one obtains the (first 
integral of the) radial geodesic equation Eq. (1!) with 
V = (1 + E 2 g tt + L 2 gM - 2 ELg^)/g rr . One can simi¬ 
larly show, by employing the definitions of the “raised” 
components (e.g., p* = mdt/dr for r proper time), that 
d<j)/dt = (g^L - g t< t‘E)/(—g tt E + g^L) and dt/dr = 
—g rr dV/dE. Using these expressions in Eq. (|A|), we 
conclude that 

d \max dt/dr dL dE 1 

dt g rr dt dt dcf>/dt 


(A3) 


When the potential admits a nearly-circular orbit (angu¬ 
lar frequency ft = d<p/dt) near the local maximum, we 
have dE « fldL in the emitted radiation and therefore 
dVjnax/dt is smaller than normal. [Similar arguments ap¬ 
ply to the minimum, and prove that stable circular orbits 
evolve to stable circular orbits.] 


APPENDIX B: KERR PARAMETERS AND 
CONSTRAINTS 


In the Teukolsky-equation-based inspiral literature 0 
0 , when the orbit is bound (=does not fall into hole or 
escape to infinity) it is characterized not by physical pa¬ 
rameters ( E,L,a ) but by the location of its radial turning 
points (r±) and the remaining root of its potential (f): 


V = 


E 2 - 1 


(r+ — r) (r — r_) (r — r) 


(Bl) 


= - (e 2 - 1) - H ^ a2 (f~ 

2 {L-aE) 2 


(B2) 


To further simplify the algebra involved, one replaces r± 
by a parameterization analogous to classical mechanics 
(semi-latus rectum and eccentricity): 


r± 


P 

1 T e 


(B3) 


(with p,e both real, positive). After replacing r± by p, e 
using r± = p/{l T e), we find the following explicit forms 
for E, L in general: (specify retrograde orbits by a nega¬ 
tive a sign) 


E 

L 


L (i ~ e2 ) 

V 2p + (1 — e 2 ) r 


sign(a) 


pVr 

\j2p+ (1 — e 2 )f 


+ aE 


(B4) 

(B5) 


The parameters p , e, f, however, are not fully inde¬ 
pendent. The coefficients of V (r ) satisfy a polynomial 
equation when the coefficients are expressed in terms of 
L,E,a', the coefficients must satisfy the same polyno¬ 
mial when the coefficients are written using r± , f. Re¬ 
expressing that polynomial in terms of p, e, we find 

0 = p 2 [p (f — 2) — 4f] 2 + a 4 [p + (l — e 2 ) f] 2 
—2a 2 p {4 (l — e 2 ) f 2 + 2 p 2 (2 + f) 

+pf [8 + (f - 2) (l - e 2 )] } (B6) 

As a practical matter, we usually specify a particle by 
(p,e,a) and then solve for f and hence all other orbital 
parameters. 


1. Separatrix 


The transition from stable to unstable occurs when 
f = r_ = p/( 1 + e). We can attempt to express this 
relation in terms of the usual orbital param eter s (p,e,a). 
Substituting in the above polynomial [Eq. (B 6 )], we get 


0 = p 2 (p s — 6 - 2e ) 2 + a 4 (e - 3) 2 (1 + e ) 2 

—2a 2 (1 + e) p s [14 + 2e 2 + p s (3 — e)] (B7) 
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where p = p s (e,a ), defined by (appropriate solutions to) 
the above relation, determines the separatrix between 
stable and unstable geodesic orbits. 


set the amplitude scale. [In this table, OT list the number 
of quadrupolar gravitational wave cycles N cyc = 2 N c .] 
We can approximate their results by the two functions 


2. Minimum of potential 


At the transition from stable to unstable orbits, the 
maximum of the potential is at r_ = r = p s /( 1 + e). 
Therefore, by symbolically differentiating Eq. (Bl), we 
find r m i n = 3p s /(3 — e). Inserting this expression into 
V and noting U(r max ) = 0, we find I max = V(r max ) - 
V(r min ): 


I -I- 


l 


(3 - e)(l + e)p s 


(B8) 


log 10 /V c ,OT ~ 0.92 + O.llOx + 0.0478x 2 , a € [0.9,0.99] 

« 1.03 + 0.073a; + 0.138a; 2 ,a < 0.9 (Cl) 

and 

lo Sio (#) OT xl0 “ 22 

« 0.3156 + 0.3193a: + 0.0451a; 2 , a e [0.8,0.99] 
« 0.191 — 0.226a; — 0.448a: 2 , a < 0.8 (C2) 


where x = log 10 (l — a). 


APPENDIX C: RESULTS FROM ORI AND 
THORNE 

So that the reader can more easily compare our esti¬ 
mates against those of Ori and Thorne, we provide ap¬ 
proximations to their results. 

In Ori and Thorne’s Table II, they tabulate S/N and 
N c for 1OM 0 on 1O 6 M 0 , using fiducial distance 1 Gpc to 
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